Systems, methods and apparatus for specialized filtered back-projection reconstruction for digital tomosynthesis

ABSTRACT

Systems, methods and apparatus are provided through which in one aspect, a three-dimensional (3D) image of an object is constructed from a plurality of two-dimensional (2D) images of the object using a specialized filter. In some embodiments, the specialized filter implements a linear ramp function, a windowing function, and/or a polynomial function. In some embodiments, the 3D image is back-projected from the filtered two-dimensional images, yielding a 3D image that has improved visual distinction of overlapping anatomic structures and reduced blurring.

RELATED APPLICATION

This application claims priority under 35 U.S.C. 120 to copending U.S.application Ser. No. 10/859,423, filed Jun. 1, 2004 entitled “SYSTEMS,METHODS AND APPARATUS FOR SPECIALIZED FILTERED BACK-PROJECTIONRECONSTRUCTION FOR DIGITAL TOMOSYNTHESIS.”

FIELD OF THE INVENTION

This invention relates generally to digital imaging, and moreparticularly to reconstructing a three dimensional image using atomosynthesis device.

BACKGROUND OF THE INVENTION

In X-ray tomosynthesis, a series of low dose X-ray images are acquiredover a range of X-ray beam orientations relative to an imaged object.Digital tomosynthesis (DTS) is a limited angle imaging technique, whichallows the reconstruction of tomographic planes on the basis of theinformation contained within the images acquired during one tomographicimage acquisition. More specifically, DTS is reconstruction ofthree-dimensional (3D) images from two-dimensional (2D) projectionimages of an object.

In DTS, one back-projection technique known as “simple back-projection”or the “shift and add algorithm” is often used to reconstruct 2D imagesinto 3D images. This technique requires a relatively straightforwardimplementation and minimal computational power requirements.

A reconstruction method used in tomosynthesis is known as the algebraicreconstruction technique (ART). Another reconstruction technique used incomputed tomography (CT) imaging (i.e., filtered back-projection)utilizes projections over the full angular range (i.e., full 360.degree.image acquisition about the object to be imaged) and a fine angularspacing between projections. Within this framework, filteredback-projection is a reconstruction method that yields high qualityreconstructions with few artifacts.

BRIEF DESCRIPTION OF THE INVENTION

The above-mentioned shortcomings, disadvantages and problems areaddressed herein, which will be understood by reading and studying thefollowing specification.

In one aspect, a method includes acquiring a plurality oftwo-dimensional X-ray images of an object from a tomosynthesis system,filtering the plurality of two-dimensional X-ray images of the objectwith a filter comprising a linear ramp function and a windowingfunction, and back-projecting the plurality of filtered two-dimensionalX-ray images into at least one three-dimensional image of the object.

In another aspect, a system includes means for acquiring a plurality oftwo-dimensional X-ray images of an object from a tomosynthesis system,means for filtering the plurality of two-dimensional X-ray images of theobject with a filter comprising a linear ramp function and a windowingfunction, and means for back-projecting the plurality of filteredtwo-dimensional X-ray images into at least one three-dimensional imageof the object.

In yet another aspect, an apparatus operable to acquire a plurality oftwo-dimensional X-ray images of an object from a tomosynthesis systemincludes apparatus operable to filter the plurality of two-dimensionalX-ray images of the object with a filter comprising a linear rampfunction and a windowing function, and apparatus operable toback-project the plurality of filtered two-dimensional X-ray images intoat least one three-dimensional image of the object.

Systems, clients, servers, methods, and computer-accessible media ofvarying scope are described herein. In addition to the aspects andadvantages described in this summary, further aspects and advantageswill become apparent by reference to the drawings and by reading thedetailed description that follows.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a block diagram that provides a system level overview of amethod to construct a three-dimensional image of an object from aplurality to two-dimensional images of the object, using specializedfilter;

FIG. 2 is a flowchart of a method of generating a three-dimensional (3D)image from two-dimensional (2D) images using specialized filter,performed by an imaging system according to an embodiment;

FIG. 3 is a flowchart of a filtering method of a specialized filterperformed by an imaging system according to an embodiment;

FIG. 4 is a diagram of an unfiltered signal of an image, according to anembodiment;

FIG. 5 is a diagram of a Ram-Lak filtered signal of an image, accordingto an embodiment;

FIG. 6 is a diagram of a Shepp-Logan filtered signal of an image,according to an embodiment;

FIG. 7 is a flowchart of a method to construct a filter by invoking aRam-Lak filter performed by an imaging system according to anembodiment;

FIG. 8 is a flowchart of a method of constructing a one-dimensional (1D)Ram-Lak filter performed by an imaging system according to anembodiment;

FIG. 9 is a flowchart of a method of filtering an image to a temporaryfile using a Ram-Lak 1D filter;

FIG. 10 is a diagram of a windowing filtered signal of an image,according to an embodiment;

FIG. 11 is a diagram of a polynomial function applied to an image,according to an embodiment;

FIG. 12 is a block diagram of the hardware and operating environment inwhich different embodiments can be practiced; and

FIG. 13 is a block diagram of the hardware and operating embodimenthaving a specialized filter.

DETAILED DESCRIPTION OF THE INVENTION

In the following detailed description, reference is made to theaccompanying drawings that form a part hereof, and in which is shown byway of illustration specific embodiments which may be practiced. Theseembodiments are described in sufficient detail to enable those skilledin the art to practice the embodiments, and it is to be understood thatother embodiments may be utilized and that logical, mechanical,electrical and other changes may be made without departing from thescope of the embodiments. The following detailed description is,therefore, not to be taken in a limiting sense.

The detailed description is divided into five sections. In the firstsection, a system level overview is described. In the second section,methods of embodiments are described. In the third section, the hardwareand the operating environment in conjunction with which embodiments maybe practiced are described. In the fourth section, particularimplementations are described. Finally, in the fifth section, aconclusion of the detailed description is provided.

System Level Overview

FIG. 1 is a block diagram that provides a system level overview of amethod to construct a three-dimensional image of an object from aplurality to two-dimensional images of the object. In method 100,projections are filtered using a specialized filter. Some embodimentsoperate in a multi-processing, multi-threaded operating environment on acomputer, such as computer 1202 in FIG. 12.

Method 100 includes filtering 102 two or more two-dimensional images ofthe object. The images are filtered in the Fourier domain using a linearramp function of the two-dimensional images, a windowing function of thetwo-dimensional images, and/or a polynomial function of thetwo-dimensional images. To combine the filtered images, the output ofthe linear ramp function, the windowing function, and/or the polynomialfunction are multiplied in the frequency domain.

The thorough and comprehensive filtering of action 102 is an improvedfiltering technique for two-dimensional images. In particular, filtering102 improves visual distinction of overlapping anatomic structures inX-ray images and reduces blurring in X-ray images. Thus filtering 102improves processing of two-dimensional images into reconstructedthree-dimensional images.

Method 100 thereafter includes back-projecting 104 the filteredtwo-dimensional images into three-dimensional images. Filtering 102 incombination with back-projecting 104 improves processing oftwo-dimensional images into reconstructed three-dimensional images.

The system level overview of the operation of an embodiment has beendescribed in this section of the detailed description. Method 100provides improved distinction of overlapping anatomic structures inX-ray images and reduces blurring in X-ray images in 3D images that arereconstructed from 2D images. While the system 100 is not limited to anyparticular image or process of back-projecting, for sake of clarity asimplified image and back-projecting has been described.

Methods of an Embodiment

In the previous section, a system level overview of the operation of anembodiment was described. In this section, the particular methodsperformed by an imaging system, a server and/or a client of such anembodiment are described by reference to a series of flowcharts.Describing the methods by reference to a flowchart enables one skilledin the art to develop such programs, firmware, or hardware, includingsuch instructions to carry out the methods on suitable computerizedsystem in which the processor of the system executes the instructionsfrom computer-accessible media. Methods 100-300 and 700-900 areperformed by a program executing on, or performed by firmware orhardware that is a part of, a computer, such as computer 1202 in FIG.12.

FIG. 2 is a flowchart of a method 200 of generating a 3D image from 2Dimages using specialized filter, performed by an imaging systemaccording to an embodiment. In method 200, a plurality of 2D views ofthe object is acquired 202. Image acquisition 202 can be performed, forexample, using any one of a number of techniques (e.g., using a digitaldetector), provided the views can be made in (or converted to) digitalform.

Amorphous silicon flat panel digital X-ray detectors are a commondetection device for tomosynthesis imaging, but in general, any X-raydetector that provides a digital projection image can be implemented,including, but not limited to, charge-coupled device (CCD) arrays,digitized film screens, or other digital detectors, such as directconversion detectors. Their low electronic noise and fast-read out timesenable acquisitions with many projections at low overall patient dosecompared to competing detector technologies.

In some embodiments, the acquired plurality of 2D views of the objectare preprocessed 204 to correction the images. Preprocessing 204 mayinclude one or more of correction for geometry effects such as distanceto X-ray tube & incident angle of X-rays at detector, correction forother system effects such as gain and offset correction, bad pixelcorrection, correction for pathlength through the tissue, taking thenegative log of the image, correction for geometry distortions (e.g.R-Square) and log transformation to restore “film-like” look, and otherpreprocessing aspects readily apparent to one of ordinary skill in theart. Preprocessing 204 may also include corrections for effects due tovarying detector positions from view to view. Preprocessing 204 may alsoinclude special pre-processing of bad detector edge correction, 2D viewweighting, and padding.

After preprocessing 204, the image value at each pixel in a viewapproximately represents the average linear attenuation value of theimaged object along the ray corresponding to that pixel, for an assumedunderlying constant thickness.

In some embodiments, each of the plurality of 2D views of the object issegmented 206. In some embodiments, segmenting 206 comprises associatingto each pixel of each 2D view, an indication of whether the pixelcontains only “air” information, or tissue information. The assigningmay implement techniques such as image value histogram segmentation,edge detection, contour following, etc. In some embodiments, segmenting206 includes the use of prior shape information (e.g., using smoothnessconstraints of the skin line), etc. In some embodiments, each view issegmented in six into pixels corresponding to rays passing through theobject and pixels corresponding to rays not passing through the object.The term “ray” refers to a given part of the X-ray beam represented by aline between the focal spot of the X-ray source and the consideredpixel. The segmenting 206 can also provide a boundary curve (i.e., acurve separating the pixels corresponding to rays passing through theobject from the 1 pixels corresponding to rays not passing through theobject), which can be used in the thickness compensation 208. Othersegmenting techniques may also be used, as would be readily apparent toone of ordinary skill in the art after reading this disclosure.

By segmenting 206, pixels not containing useful information (i.e.,pixels corresponding to rays not passing through the object) can begiven a constant image value, in some embodiments, about equal to anaverage image value for pixels corresponding to rays passing through theobject. This enhances the appearance of structures (includingabnormalities) within the reconstructed 3D image of the object, andreduces artifacts. Thus, the overall performance of the reconstructioncan be greatly improved. As will be described in detail below,segmenting 206 the 2D views is a particularly effective technique whenused in combination with the thickness compensation step 208.

Method 200 also includes compensating 208 the segmented 2D views of theobject for a thickness of the image object. Conventional thicknesscompensation techniques can be used in the segmenting 208. Compensating208 for thickness provides a significant reduction of reconstructionartifacts due to the reduced thickness of the imaged object near theboundary and preserves coarse scale in the image corresponding tovariations in tissue characteristics that are not due to the reducedthickness.

Compensating 208 for thickness allows for a “fair comparison” ofdifferent image values in the back-projected 2D views, in thatback-projecting 104 in some embodiments uses an order statisticsoperator, and therefore compares different image values from differentprojection images. A bias in one or more values which is due to areduced thickness at the corresponding location can have acounterproductive effect on the resulting reconstructed 3D image. Thiseffect can be minimized by compensating 208 for thickness. Thus,compensating 208 for thickness provides substantial improvements overconventional techniques.

The plurality of 2D views are filtered 102. Filtering 102 can preferablybe implemented as a one-dimensional (1D) or a 2D filtering process. In1D filtering, appropriate filters are mapped from CT geometry to atomosynthesis geometry. One example of CT geometry is 360 degreeacquisition, in which the detector rotates opposite of the tube suchthat the incoming X-rays are substantially perpendicular to the detectorsurface for all views. One example of tomosynthesis geometry is lessthan 360 degree acquisition, and the angle of the incoming X-rays on thedetector varies from view to view. In some embodiments, the mappedfilter would be shift-variant, but the effect on the image quality ofthe reconstructed volume is generally negligible. Mapped versions ofconventional filters, as well as other 1D filters which are optimizedwith respect to some reconstruction image quality criterion can beimplemented. Filters for each view can be used in accordance withspecific acquisition geometry. The 2D filter can be generated from theone-dimensional filter by either swirl it 360 degrees (“circular”), ormultiplying itself by its transposed version (“rectangular”).

In some embodiments of method 200, constraints are applied 210 to thefiltered 2D views. In more specific embodiments of applying 210constraints, only “physically admissible” image values are retained forfurther analysis. For example, negative values (which do not correspondto physical reality) may be set to zero, or the maximum attenuation ofthe material of the imaged object may be known, which would allow one toderive a maximum meaningful value, and the image could thus be truncatedto that maximum meaningful value. As only physically admissible imagevalues are retained, constraint application improves the noise andartifact characteristics of the reconstruction method.

The filtered plurality of 2D views of the object are then back-projected104 into a 3D representation of the object. In some embodiments theback-projecting 104 uses an order statistics-based back-projectingtechnique.

Order statistics-based back-projection is significantly different inmany aspects to simple back-projection reconstruction. Specifically, inorder statistics based back-projecting, the averaging operator which isused to combine individual back-projected image values at any givenlocation in the reconstructed volume is replaced by an order statisticsoperator. Thus, instead of simply averaging the back-projected pixelimage values at each considered point in the reconstructed volume, anorder statistics based operator is applied on a voxel-by-voxel basis.

Depending on the specific framework, different order statisticsoperators may be used (e.g., minimum, maximum, median, etc.), but insome embodiments, an operator which averages all values with theexception of some maximum and some minimum values is preferred. Moregenerally, an operator which computes a weighted average of the sortedvalues can be used, where the weights depend on the ranking of theback-projected image values. In particular, the weights corresponding tosome maximum and some minimum values may be set to zero. By using theaforementioned operator for breast imaging, streak artifacts (which aregenerally caused either by high contrast structures—maxima, or by the“overshoot” caused by the filtering of some high contraststructure—minima) are minimized, while some of the noise reductionproperties of the linear averaging operator are retained.

Alternatively, other back-projection methods may also be used in action104, such as Shift & Add, Generalized Filtered Back-Projection. Othermethods may also be included, based on a minimum-norm solution such asART, DART, MITS, TACT, Fourier-Based Reconstruction, ObjectiveFunction-Based Reconstruction, or combinations thereof.

Back-projecting 104 is further improved by back-projecting data alreadysegmented in action 206. For example, the segmentation result can beused to set a reconstructed value to zero (or some other appropriatevalue) if at least a single back-projected image value indicates an“outside” location (i.e., the corresponding pixel in that view wasdetermined not to correspond to a ray passing through the imagedobject). In addition, if some voxel in the reconstruction volume is notcontained in all projection radiographs (e.g., because for someprojection angle the corresponding point was not projected onto theactive area of the detector), then only the projection radiographs thatcontain this voxel are taken into account in the reconstruction. In analternate approach, one can artificially increase the image size byadding regions to the boundaries of the image and setting the imagevalues in these regions equal to the “background value.” Both of theseapproaches help to minimize artifacts which are due, for example, to theboundary effects of the detector.

In some embodiments, method 200 includes applying 212 constraints to thereconstructed dataset/3D image after the back-projecting 104. Applyingconstraints 212 may include setting negative values to zero, truncatinghigh values to the maximum value for the type of object imaged, etc.Apply 212 constraints may also include post-processing functions such asimage/contrast enhancement of tissue equalization, thicknesscompensation, and/or brightness/white balancing. This may beparticularly useful in combination with iterative update step 216described in detail below.

Thereafter, the reconstructed 3D image of the object is output 214, suchas by displaying the image.

In some embodiments, the reconstructed 3D representation of the objectis iteratively updated 216. In various embodiments, iteratively updating216 includes re-projecting the reconstructed 3D image of the object,comparing the reconstructed 3D image to the acquired views, and updatingthe reconstructed 3D image of the object. In some embodiments,iteratively updating 216 is performed prior to outputting 214 thereconstructed 3D image of the object. In some embodiments, iterativelyupdating 216 is performed at intervals followed by an updated outputted3D image.

FIG. 3 is a flowchart of a filtering method 300 of a specialized filterperformed by an imaging system according to an embodiment. Method 300 isone embodiment of filtering 102 in FIG. 1.

In method 300, a linear function 302, a windowing function 304 and apolynomial function are performed on the 2D images. The combination ofthe three filtering actions 302, 304, and 306 provides a thorough andcomprehensive filtering that improves distinction of overlappinganatomic structures and reduces blurring in 3D images that arereconstructed from 2D images that have been filtered in accordance withmethod 300.

The linear function 302 provides high-pass filtering that deemphasizeshigh frequencies in which the gain is proportionate to the frequency. Insome embodiments, the linear function is a starter kernel function. Insome embodiments, the linear function is a Ramachandran-Lakshminarayanan(Ram-Lak) filter or a Shepp-Logan filter. A Ram-Lak filter is also knownas a ramp function. The Ram-Lak filter function shown in Formula 1below:H(ξ)=|ξ|rect(ξ/(2ξ_(max)))

Formula 1

In Formula 1, ξ is frequency. Frequency components greater than ξ_(max)are truncated. A conventional example of a filter that reduces noise bysuppressing the gain in high frequencies is the Shepp-Logan filter,shown in Formula 2 below:H(ξ)=|ξ|sinc(ξ/(2ξ_(max)))rect(ξ(2ξ_(max)))

Formula 2

The Shepp-Logan filter shown in Formula 2 is a modified version of theRam-Lak filter. The Shepp-Logan filter multiplies the Ram-Lak filter bythe sinc function, which is the equivalent of the convolution with therect function in the real domain, which has the effect of averagefiltering in the real domain.

The windowing function 304 can be a Hanning function or any otherbell-shaped Gaussian function that ramps down smoothly.

FIG. 4 is a diagram 400 of an unfiltered signal 402 of an image,according to an embodiment. The diagram shows the unfiltered signal 402having a frequency ξ 404 plotted along an amplitude of H(ξ) 406.

FIG. 5 is a diagram 500 of a Ram-Lak filtered signal 504 of an image,according to an embodiment. The diagram shows the signal 502 that isfiltered by a Ram-Lak filter in accordance with Formula 1 having afrequency ξ 504 plotted along an amplitude of H(ξ) 506. The portion 508of the signal 504 that is outside the bounds −ξ_(max) 510 and ξ_(max)512 is truncated.

FIG. 6 is a diagram 600 of a Shepp-Logan filtered signal 602 of animage, according to an embodiment. The diagram shows the signal 602 thatis filtered by a Shepp-Logan filter in accordance with Formula 2 havinga frequency ξ 604 plotted along an amplitude of H(ξ) 606. The portion608 of the signal 602 that is outside the bounds −ξ_(max) 610 andξ_(max) 612 is truncated. The slope k 514 is optimized and balanced toimprove contrast and minimize image noise.

FIG. 7 is a flowchart of a method 700 to construct a filter by invokinga Ram-Lak filter performed by an imaging system according to anembodiment. Method 700 is one embodiment of a liner filtering function302 in FIG. 3.

Method 700 includes initializing 702 input and output files. In someembodiments, that includes creating an output filename, opening andinput image file for reading, opening the output file for writing,reading in a header, obtaining the number of images and image size fromthe header, and writing the header to the output file.

After initialization, a 1D filter is constructed 704 in the fast-Fouriertransform (FFT) domain, such as by invoking Ram-Lak 1D filterconstructor method 800 in FIG. 8. In other embodiments, other filtersare constructed. The 1D filter implements a data structure thatrepresents a vector having 4096 elements. The size of 4096 elements isdetermined using the Nyquist criteria. In 1927, Mr. Nyquist of Bell Labsdetermined that an analog signal should be sampled at least twice thefrequency of its highest-frequency component in order to be convertedinto an adequate representation of the signal in digital form. This ruleis now known as the Nyquist-Shannon sampling theorem. In digitaltomosynthesis, a detector will generate a digital representation of animage using 2022 elements. Using a Nyquist criteria of a multiple of 2,a minimum of 4044 elements is required to adequately represent theimage. The number of 4044 elements is under up to 4096 in order to bereadily addressed in binary addressing that is used by computers.

Thereafter, the image is filtered 706 to a temporary file having 4096elements using the constructed Ram-Lak 1D filter from action 704. Oneexample of the filtering 706 is method 900 in FIG. 9.

Subsequently, method 700 includes receiving 708 images from thetemporary file and applying 710 a scaling factor to the temporary file.The temporary file is then written 712 to the output file.

FIG. 8 is a flowchart of a method 800 of constructing a 1D Ram-Lakfilter performed by an imaging system according to an embodiment. Method800 is one embodiment of action 704 in FIG. 7.

Method 800 includes determining 802 a polynomial tweaking windowcoefficients and a cut-off frequency. Method 800 also includesdetermining 804 a start-array index and an end-array index for windowingramp function. Method 800 also includes filling 806 a basic window withbinary ‘1’'s up to and including an element following the start-arrayindex. Method 800 subsequently includes multiplying 808 a polynomialwindow. Thereafter, method 800 includes copying 810 the left half windowto the right half of the window.

FIG. 9 is a flowchart of a method 900 of filtering 706 an image to atemporary file using a Ram-Lak 1D filter. Method 900 is one embodimentof filtering 706 in FIG. 7.

Method 900 includes determining 902 a correct filter scaling factor.Then a determination 904 is made as to whether or not additional imagesare to be processed.

If more images are to be processed, then the next image is placed 906 inthe FFT domain. Then a filter is applied 908 to that image, and a bandof values of the image are identified 910. The image is saved 912 in atemporary file.

When no more images are to be processed, the temporary file is closed914.

FIG. 10 is a diagram 1000 of a windowing filtered signal 1002 of animage, according to an embodiment. Diagram 1000 shows the signal 1002that is filtered by a windowing function filter in accordance withaction 304 in FIG. 3 plotted along a magnitude 1004 and a frequency1006. The windowing function is a continuous function that has a plateausegment 1008 before it starts to ramp down smoothly from a maximummagnitude of 1.0 1010 to a magnitude of 0.0 (zero) 1012 at a cutofffrequency (fc) 1014. In other embodiments, the function starts at amagnitude 1004 value of other than 1.0 and does not have a plateausegment before it starts to ramp down. A starting frequency (fs) 1016identifies the beginning of the plateau segment. In some embodiments,the cutoff frequency (fc) 1014 and the starting frequency (fs) 1016 areoptimized to preserve image details and produce minimal image noises.

FIG. 11 is a diagram 1100 of a polynomial function 1002 applied to animage, according to an embodiment. Diagram 1100 shows the polynomialfunction 1102 in accordance with action 306 in FIG. 3 plotted along amagnitude 1104 and a frequency 1106. The polynomial function 1100 finetunes the frequency response of the specialized filter to enhancecertain frequencies and to suppress others. The polynomial function istypically a continuous function of second-order (not shown) orhigher-order polynomials, such as the 4^(th)-order polynomial shown indiagram 1100. An example of a fourth-order polynomials is shown inFormula 3 below:F(w)=c4w ⁴ +c3w ³ +c2w ² +c1w+c ₀

Formula 3

The fourth-order polynomial of Formula 3 is shown in dotted-lines inFIG. 11, and the product of the polynomial functions and a windowingfunction are also shown in solid-lines. Another term for a polynomialfunction includes a piece-wise spline function. Applying a polynomialfunction helps to differentiate the particular anatomies/structure inthe image.

In some embodiments, methods 100-300 and 700-900 are implemented as acomputer data signal embodied in a carrier wave, that represents asequence of instructions which, when executed by a processor, such asprocessor 1204 in FIG. 12, cause the processor to perform the respectivemethod. In other embodiments, methods 100-300 and 700-900 areimplemented as a computer-accessible medium having executableinstructions capable of directing a processor, such as processor 1204 inFIG. 12, to perform the respective method. In varying embodiments, themedium is a magnetic medium, an electronic medium, or an optical medium.

Hardware and Operating Environment

FIG. 12 is a block diagram of the hardware and operating environment1200 in which different embodiments can be practiced. The description ofFIG. 12 provides an overview of computer hardware and a suitablecomputing environment in conjunction with which some embodiments can beimplemented. Embodiments are described in terms of a computer executingcomputer-executable instructions. However, some embodiments can beimplemented entirely in computer hardware in which thecomputer-executable instructions are implemented in read-only memory.Some embodiments can also be implemented in client/server computingenvironments where remote devices that perform tasks are linked througha communications network. Program modules can be located in both localand remote memory storage devices in a distributed computingenvironment.

Computer 1202 includes a processor 1204, commercially available fromIntel, Motorola, Cyrix and others. Computer 1202 also includesrandom-access memory (RAM) 1206, read-only memory (ROM) 1208, and one ormore mass storage devices 1210, and a system bus 1212, that operativelycouples various system components to the processing unit 1204. Thememory 1206, 1208, and mass storage devices, 1210, are types ofcomputer-accessible media. Mass storage devices 1210 are morespecifically types of nonvolatile computer-accessible media and caninclude one or more hard disk drives, floppy disk drives, optical diskdrives, and tape cartridge drives. The processor 1204 executes computerprograms stored on the computer-accessible media.

Computer 1202 can be communicatively connected to the Internet 1214 viaa communication device 1216. Internet 1214 connectivity is well knownwithin the art. In one embodiment, a communication device 1216 is amodem that responds to communication drivers to connect to the Internetvia what is known in the art as a “dial-up connection.” In anotherembodiment, a communication device 1216 is an Ethernet® or similarhardware network card connected to a local-area network (LAN) thatitself is connected to the Internet via what is known in the art as a“direct connection” (e.g., T1 line, etc.).

A user enters commands and information into the computer 1202 throughinput devices such as a keyboard 1218 or a pointing device 1220. Thekeyboard 1218 permits entry of textual information into computer 1202,as known within the art, and embodiments are not limited to anyparticular type of keyboard. Pointing device 1220 permits the control ofthe screen pointer provided by a graphical user interface (GUI) ofoperating systems such as versions of Microsoft Windows®. Embodimentsare not limited to any particular pointing device 1220. Such pointingdevices include mice, touch pads, trackballs, remote controls and pointsticks. Other input devices (not shown) can include a microphone,joystick, game pad, satellite dish, scanner, or the like.

In some embodiments, computer 1202 is operatively coupled to a displaydevice 1222. Display device 1222 is connected to the system bus 1212.Display device 1222 permits the display of information, includingcomputer, video and other information, for viewing by a user of thecomputer. Embodiments are not limited to any particular display device1222. Such display devices include cathode ray tube (CRT) displays(monitors), as well as flat panel displays such as liquid crystaldisplays (LCD's). In addition to a monitor, computers typically includeother peripheral input/output devices such as printers (not shown).Speakers 1224 and 1226 provide audio output of signals. Speakers 1224and 1226 are also connected to the system bus 1212.

Computer 1202 also includes an operating system (not shown) that isstored on the computer-accessible media RAM 1206, ROM 1208, and massstorage device 1210, and is and executed by the processor 1204. Examplesof operating systems include Microsoft Windows®, Apple MacOS®, Linux®,UNIX®. Examples are not limited to any particular operating system,however, and the construction and use of such operating systems are wellknown within the art.

Embodiments of computer 1202 are not limited to any type of computer1202. In varying embodiments, computer 1202 comprises a PC-compatiblecomputer, a MacOS®-compatible computer, a Linux®-compatible computer, ora UNIX®-compatible computer. The construction and operation of suchcomputers are well known within the art.

Computer 1202 can be operated using at least one operating system toprovide a graphical user interface (GUI) including a user-controllablepointer. Computer 1202 can have at least one web browser applicationprogram executing within at least one operating system, to permit usersof computer 1202 to access intranet or Internet world-wide-web pages asaddressed by Universal Resource Locator (URL) addresses. Examples ofbrowser application programs include Netscape Navigator® and MicrosoftInternet Explorer®.

The computer 1202 can operate in a networked environment using logicalconnections to one or more remote computers, such as remote computer1228. These logical connections are achieved by a communication devicecoupled to, or a part of, the computer 1202. Embodiments are not limitedto a particular type of communications device. The remote computer 1228can be another computer, a server, a router, a network PC, a client, apeer device or other common network node. The logical connectionsdepicted in FIG. 12 include a local-area network (LAN) 1230 and awide-area network (WAN) 1232. Such networking environments arecommonplace in offices, enterprise-wide computer networks, intranets andthe Internet.

When used in a LAN-networking environment, the computer 1202 and remotecomputer 1228 are connected to the local network 1230 through networkinterfaces or adapters 1234, which is one type of communications device1216. Remote computer 1228 also includes a network device 1236. Whenused in a conventional WAN-networking environment, the computer 1202 andremote computer 1228 communicate with a WAN 1232 through modems (notshown). The modem, which can be internal or external, is connected tothe system bus 1212. In a networked environment, program modulesdepicted relative to the computer 1202, or portions thereof, can bestored in the remote computer 1228.

Computer 1202 also includes power supply 1238. Each power supply can bea battery.

FIG. 13 is a block diagram of the hardware and operating embodiment 1300having a specialized filter. Computer 1202 includes a linear functioncomponent 1302 that implements the linear function of action 302 in FIG.3. Computer 1202 also includes a windowing function component 1304 thatimplements the windowing function of action 304 in FIG. 3. Computer 1202further includes a polynomial function component 1306 that implementsthe polynomial function of action 306 in FIG. 3.

Components 1240, 1242 and 1244 and components that implement methods100-300 and 700-900 can be embodied as computer hardware circuitry or asa computer-accessible program, or a combination of both. In anotherembodiment, the components are implemented in an application serviceprovider (ASP) system.

More specifically, in the computer-accessible program embodiment, theprograms can be structured in an object-orientation using anobject-oriented language such as Java, Smalltalk or C++, and theprograms can be structured in a procedural-orientation using aprocedural language such as COBOL or C. The software componentscommunicate in any of a number of means that are well-known to thoseskilled in the art, such as application program interfaces (API) orinterprocess communication techniques such as remote procedure call(RPC), common object request broker architecture (CORBA), ComponentObject Model (COM), Distributed Component Object Model (DCOM),Distributed System Object Model (DSOM) and Remote Method Invocation(RMI). The components execute on as few as one computer as in computer1202 in FIG. 12, or on at least as many computers as there arecomponents.

CONCLUSION

Systems, methods and apparatus that generate a 3D image from 2Dtomosynthesis images using a specialized filter has been described.Although specific embodiments have been illustrated and describedherein, it will be appreciated by those of ordinary skill in the artthat any arrangement which is calculated to achieve the same purpose maybe substituted for the specific embodiments shown. This application isintended to cover any adaptations or variations. For example, althoughdescribed in procedural design terms, one of ordinary skill in the artwill appreciate that implementations can be made in an object-orienteddesign environment or any other design environment that provides therequired relationships.

In particular, one of skill in the art will readily appreciate that thenames of the methods and apparatus are not intended to limitembodiments. Furthermore, additional methods and apparatus can be addedto the components, functions can be rearranged among the components, andnew components to correspond to future enhancements and physical devicesused in embodiments can be introduced without departing from the scopeof embodiments. One of skill in the art will readily recognize thatembodiments are applicable to future communication devices, differentfile systems, and new data types.

The terminology used in this application is meant to include all imageand communication environments and alternate technologies which providethe same functionality as described herein.

1. A method comprising: acquiring a plurality of two-dimensional X-rayimages of an object from a tomosynthesis system; filtering the pluralityof two-dimensional X-ray images of the object with a filter comprising alinear ramp function and a windowing function; back-projecting theplurality of filtered two-dimensional X-ray images into at least onereconstructed three-dimensional image of the object; and iterativelyupdating the reconstructed three dimensional image of the object.
 2. Themethod of claim 1, wherein the filtering further comprises filteringperformed in one dimension.
 3. The method of claim 1, wherein thefiltering further comprises filtering performed in two dimensions. 4.The method of claim 1, wherein the linear ramp function is performedbefore the windowing function and a polynomial function.
 5. The methodof claim 1, wherein the linear ramp function is performed before apolynomial function.
 6. The method of claim 1, wherein the linear rampfunction further comprises a Ram-Lak filter.
 7. The method of claim 1,wherein the linear ramp function further comprises a Shepp-Logan filter.8. The method of claim 1, wherein the windowing function farthercomprises a Hanning filter.
 9. The method of claim 1 further comprising:before the filtering: preprocessing the plurality of the two-dimensionalimages of the object; segmenting each of the plurality of thepreprocessed two-dimensional images of the object; and compensating thesegmented two-dimensional images of the object.
 10. The method of claim1 further comprising: after the back-projecting: outputting thereconstructed three dimensional image of the object.
 11. A systemcomprising: apparatus operable to acquire a plurality of two-dimensionalX-ray images of an object from a tomosynthesis system; apparatusoperable to filter the plurality of two-dimensional X-ray images of theobject with a filter comprising a linear ramp function and a windowingfunction; apparatus operable to back-project the plurality of filteredtwo-dimensional X-ray images into at least one reconstructedthree-dimensional image of the object; and apparatus operable toiteratively update the at least one reconstructed three-dimensionalimage of the object.
 12. The system of claim 11 further comprising:apparatus operable to preprocess the plurality of the two-dimensionalX-ray images of the object; apparatus operable to segment each of theplurality of the preprocessed two-dimensional X-ray images of theobject; and apparatus operable to compensate the segmentedtwo-dimensional X-ray images of the object.
 13. The system of claim 11further comprising: apparatus operable to apply constraints to thefiltered two-dimensional images of the object.
 14. The system of claim11 further comprising: apparatus operable to filter the plurality oftwo-dimensional X-ray images of the object from a polynomial function ofthe two-dimensional X-ray images; and apparatus operable to output thereconstructed three-dimensional image of the object.
 15. The system ofclaim 11 further comprising: apparatus operable to digital detection ofX-ray radiation.
 16. The system of claim 11 further comprising:apparatus operable to preprocess the plurality of the two-dimensionalX-ray images of the object; apparatus operable to segment each of theplurality of the preprocessed two-dimensional X-ray images of theobject; apparatus operable to compensate the segmented two-dimensionalX-ray images of the object; apparatus operable to apply constraints tothe filtered two-dimensional images of the object; apparatus operable tofilter the plurality of two-dimensional X-ray images of the object froma polynomial function of the two-dimensional X-ray images; apparatusoperable to output the reconstructed three-dimensional image of theobject; and apparatus operable to iteratively update the reconstructedthree-dimensional image of the object.
 17. A device comprising:apparatus operable to acquire a plurality of two-dimensional X-rayimages of an object from a tomosynthesis system; apparatus operable tofilter the plurality of two-dimensional X-ray images of the object witha filter comprising a linear ramp function and a windowing functionapparatus operable to back-project the plurality of filteredtwo-dimensional X-ray images into at least one reconstructedthree-dimensional image of the object; apparatus operable to iterativelyupdate the at least on reconstructed three-dimensional image of theobject.
 18. The device of claim 17 further comprising: apparatusoperable to preprocess the plurality of the two-dimensional X-ray imagesof the object; apparatus operable to segment each of the plurality ofthe two-dimensional X-ray images of the object; and apparatus operableto compensate the segmented two-dimensional X-ray images of the object.19. The device of claim 17 further comprising: apparatus operable toapply constraints to the filtered two-dimensional images of the object.20. The device of claim 17 further comprising: apparatus operable tofilter the plurality of two-dimensional X-ray images of the object froma polynomial function of the two-dimensional X-ray images; and apparatusoperable to display the at least one reconstructed three-dimensionalimage of the object.